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Abstract 

We use a novel approach to analyze the one-dimensional spinless Falicov-Kimball model. We 
derive a simple effective model for the occupation of the localized orbitals which clearly reveals 
the origin of the known ordering. Our study is extended to a quantum model with hybridization 
between the localized and itinerant states: we find a crossover between the well-known weak- 
and strong-coupling behaviour. The existence of electronic polarons at intermediate coupling is 
confirmed. A phase diagram is presented and discussed in detail. 

PACS numbers: 71.10.Fd, 71.28.+d, 71.27.+a, 71.30.+h 
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The Falicov-Kimball model (FKM) is one of the simplest yet most versatile models of 
strongly correlated electron systems. It describes a band of conduction (o) electrons inter- 
acting via a repulsive contact potential G with localized /-electrons of energy tf. Originally 
conceived of as a model for valence transitions (VTs) Q|, the FKM is today usually regarded 
as a simple model of a binary alloy, the so-called crystallization problem (CP) [2]. The CP 
assumes fixed c- and /-electron populations, examining the diagonal LRO (or "crystalline 
order") adopted by the /-electrons in the ground state. In spite of detailed numerical maps 
of the one- dimensional (ID) FKM's phase diagram Q], its behaviour is only rigorously un- 
derstood at weak- jj| and strong-coupling 0. 

A c-f hybridization term is usually added to the FKM in order to model mixed-valence 
systems j6|. The presence of quantum valence fluctuations makes this a fundamentally 
different problem to the FKM: it is referred to as the extended or quantum Falicov-Kimball 
Model (QFKM). In contrast to the FKM, the QFKM's behaviour remains largely unknown. 

nn 

Although there is no consensus on the nature of the ground state pj,l2|, excitonic effects are 
expected to dominate the physics and so this remains an important open question in the 
theory of optical properties of strongly correlated electron systems. 

Using the well-known non-perturbative technique of bosonization, we obtain for the first 
time a complete theoretical description of the ID FKM below half-filling. This analysis 
reveals the role of the c-electrons in generating the observed /-electron order; all aspects of 
the CP results are naturally explained within our approach. We derive an effective Hamilto- 
nian for the /-orbital occupancy which accurately predicts the FKM's phase diagram. For 
a small hybridization potential our approach is also suitable for the QFKM. We present a 
phase diagram that interpolates between the well-known weak- and strong-coupling limits; 
our work rigorously establishes the relevance of electronic polaron effects. 

The FKM for spinless fermions has the Hamiltonian 

w FKM = -t c \ c i + e f^2 n j + G J2 nC j n j W 
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We consider only repulsive potentials, G > 0. The concentration of electrons is fixed at 
n = jj: J2j{ n j + n< j} < 1 where N is the number of sites. In our study of the QFKM, we 
adopt the standard on-site hybridization potential H hyb = VJ2j{fj c j + H.c.}. 

Proceeding with our bosonization solution, we linearize the c-electron spectrum about 
the two Fermi points and define left- and right-moving Fermion fields c u j, v = L{— ), R(+) 



respectively as subscript (otherwise). The chiral density operators p u (k) = J2w c lk'+k c ^k' 
obey the standard Luttinger commutators [p v {k), p v >{k')]_ = 8 u y5k-k'vkL/2ir for a system 
of size _L 3> a, the lattice constant. The dual Bose fields are constructed in terms of the 
Pu{k) 



= »T£^M*)A«(*)e <to ' (3) 
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The cut-off function A a {k) appearing in Eq. (J2J and Eq. Q satisfies the conditions A a (k) ~ 1 
for \k\ < - and A a (k) « otherwise. A Q (k) enforces the finite minimum wavelength a > a of 
the bosonic density fluctuations p v {k): taking into account this wavelength limit is essential 
for preserving the lattice structure of the FKM Q|. Since the Bose fields cannot 'resolve' 
distances less than a, we find that the usual field commutators are 'smeared' over this length 
scale: 

[(j)(Xj),e(x f )]_ = —Sgll^Xy-Xj) (4) 

[d x (j)(x j ),9(x j i)]_ = -in5 a (xji - Xj) (5) 



Here sgn. a {x) and S a (x) are the a-smeared sign and Dirac delta functions respectively 11 1. 
The chiral Fermions are represented in terms of the Bose fields by the Mandelstram iden- 



tity c V j = y Aa/aF u exp (—iv [(J)(xj) — v9(xj)\). The dimensionless normalization constant 
A is dependent upon the form of A a (/c); the Klein factors F v are "ladder operators" between 
subspaces of differing c-electron number. The bosonic form of the number operators is given 
by 

7r a „ 

The first term on the RHS, Uq, is the non-interacting c-electron concentration; the second 
term is the forward-scattering density fluctuation; the third term is the first order backscat- 
tering correction. Higher order corrections are usually neglected. 

Our bosonization analysis requires that Uq ^ 0. For the sake of brevity, here we as- 
sume the symmetric case with equal c- and /-populations in the non-interacting limit, i.e. 
n o = n o = \ n - F° r the VT problem, this pins the Fermi level at e/ = — 2t cos(kpa), 
kF = irn/2a is the Fermi momentum (note that e/ is irrelevant to the CP). Using the 
electron concentration condition we also re- write the Coulomb interaction up to a constant: 



G J2j n< j n j — G J2j n j (jij ~ |) — \G J2j n j- K is now a simple matter to obtain the bosonized 
form of the Hamiltonian by substituting these identities into Eq. (JTJ. 

We apply a shift transformation on the c-electron bosonic fields U = 
exp {i-^: J2j{ n j ~ ^)9{ x j)\ wnere vp = 2tasm(k F a) is the c-electron Fermi velocity. 
This transformation rotates the basis of the Hilbert space so that the c-electrons are 
explicitly coupled to the /-orbitals. Introducing a pseudospin-i representation for the 
/-electron occupation at site j, nj — | = rj, we write the transformed Hamiltonian 
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t* cos ^2 ^>(xj) - K a {xj) - k F + ^- xXj (7) 



where K. a {xj) = S a {xj) + C a (xj), S a {xj) = f E^=i[sgn a (x i+n - Xj) - 1] (r/ +n - r/_ n+1 ) 
and C a {xj) = § - l) £~ =1 sgn Q (x i+n - x-) (r? +n - r/_ n J. S a (xj) is a measure of the 
short-range behaviour of the pseudospins, subtracting the t z within a to the left of xj from 
the t z within a to the right; C a (xj) probes the long-range behaviour, subtracting the r z 



12 1 . Note that we retain 



more than a to the left of Xj from the r z more than a to the right 
all terms produced by the canonical transform in Eq. (J7|). 

The removal of the forward-scattering Coulomb interaction introduces the important new 
term 

G 2 a 2 

n int = -- — T ! 6 *(xj - Xj'Y], (8) 

2rrv F . ., 

The physical origin of this interaction is the derealization of the c-electrons below scales 
less than a > a. A c-electron spread over more than a single lattice site carries the same 
charge over these sites: due to the Coulomb interaction G, this will favour empty underlying 
/-orbitals. Although such "segregation" is known to occur in the CP 

the derivation 

of the responsible effective interaction is a new achievement in the history of the model. 
The a-smeared 5-function in Eq. (jHJ) makes it clear that the interaction is short-ranged; we 
may therefore approximate the interaction as nearest-neighbour only, H int ~ —J'J2jT?Tj +1 , 
J = G 2 a 2 5 a (a)/7cv F . 

The other terms in Eq. Q involving the pseudospins are a constant and a site-dependent 
longitudinal field. The former is important only in the VT: together with the first term 
in Eq. (jJJ) this determines the distribution of the total electron population across the two 



4 



orbitals. The site-dependent field arises from the 2/cp-backscattering of the c-electrons off 
the /-orbitals. This is responsible for the observed LRO phases, in close analogy to the 
Pe.eHs state Q. It , clear from Eq. Q that the LRO donates the segregation at weak- 
coupling G < f; increasing G, however, Eq. (jHJ) eventually causes the system to segregate 
(the SEG phase). 

The competition between the SEG and LRO phases can be simply studied within the 
framework of our bosonization approach: since the only coupling in Eq. (J2J) between the Bose 
fields and the r-pseudospins is in the backscattering term, by replacing <fi with a suitably 
chosen expectation value we obtain an effective Hamiltonian 7i cS for the /-occupation only. 
Exact diagonalization calculations on 3200-site chains reveals that within the SEG phase 



the c-electrons are at their Luttinger liquid fixed point for all G 



is therefore valid across the phase diagram |10j . Substituting this into Eq. (JJJ), we find 



13]. 



The choice (<f>(xj)) = 



Ti.eS — J X] T j T j + 1 
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r* cos (2 {K a (xj) + [k F + Xj }) 

j 

(9) 

Below we present a detailed analysis of Ti, eS . 

Since the c- and /-populations are fixed in the CP, Eq. (jHJ) by itself fully describes the 
system close to the SEG-LRO boundary. To determine this line the site-dependence of the 
longitudinal field must be known; in particular, we examine SrJxj) and C a (xj). Within the 
SEG phase, the /-electrons are arranged into a single block [a]. The magnitude of C a (xj) 
reaches a maximum at the edge of this block, increasingly linearly as the edge is approached 
from either side. In contrast, S a (xj) vanishes everywhere except in the vicinity of this edge. 
The longitudinal field thus has the general form hj oc cos(uj + </>), u and constants. 

The periodicity (u) of M determines the critical Ising coupling sl4]: for G t we find 



the critical line G c a = vf\J2tiA/ a5 a {a). The /-electron configuration acts as an applied 
scattering potential G(tj) for the c-electrons ^j]. Below the average interparticle separation 
~ kp 1 the c-electrons scatter independently: for n<l the low-energy (E = 0) wavefunctions 
if) therefore obey the Schrodinger equation d^ifj(x) = Gm(r^)i/j(x) in the continuum limit 
(m is the bare electron mass). From elementary quantum mechanics we find plane wave 
solutions over unoccupied orbitals ((t*) = — |) and exponentially decaying solutions with 
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characteristic length ( oc yt/G over occupied orbitals ((t*) = |) [15j . ( clearly corresponds 
to a in the boson theory; for small G c Taylor-expanding 5 a (a) in the critical line equation 
gives + a 2 /2a 2 + . . .), implying a linear relationship between G c a/vF 

and G c . This is confirmed by examining the numerical data of Gajek et al, with our fit 
giving A = 1/327T 2 , a = a/2y/GZ, and a critical line G c /t = 0.5 sin(mr/2)/(l — sin(n7r/2)) 
as plotted in Fig. This is the first analytic expression for the SEG-LRO boundary valid 
at weak coupling. 

Moving on to the VT problem, it is clear from Eq. J7J) that the Coulomb interaction 
will shift the /-level downwards, emptying the c-band. The depletion of the c-electron 
population causes segregation to occur at a smaller value of G than in the CP- the new 
boundary can be calculated self-consistently from the available numerical results [3J]. Using 
the same curve-fitting technique as for the CP, we find excellent agreement with the data 
for GJt = 0.35sin(wr/2)/[l - sin(mr/2)] (A « 1.55 x 10~ 3 , a « 0.418a/ see Fig. flj. 
Further increasing G above the critical value G/t = 4(1 — cos(mr/2))/(l — n) the /-level 
lies below the c-electron band edge and so all electrons have /-character. The absence of 
c-electrons to mediate the segregating interaction implies that here any /-configuration is 
the ground state. 

We extend our analysis to the physically interesting limit V<St of the QFKM. Bosonizing 
and carrying out the canonical transform, the hybridization adds the term 



to the effective pseudospin Hamiltonian Eq. (|5|). We have used a generalized Jordan- Wigner 
transformation to combine the /-electron operators and the Klein factors into pseudospins 



dospin Hamiltonian is easily recognized as the Ising model in site-dependent longitudinal 



Although we are not aware of any systematic study of this Hamiltonian, we can never- 
theless deduce important aspects of its physical behaviour. The two extreme limits G~^>t 
and G<^Lt are easiest to discuss. In the former case the presence of V <C G does not rad- 
ically modify the VT physics. Where the longitudinal field dominates, we can discard the 
Ising interaction so that the pseudospins track the longitudinal and transverse fields; clearly 
the transverse field will only be important where the longitudinal field is small. This will 




(10) 




The resulting effective pseu- 



(hj) and transverse (/ij) fields. 
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"smear" the crystalline arrangement, but leaves the basic picture of a density wave intact. 
Where the Ising interaction dominates, the hybridization will stabilize the SEG phase as 
quantum fluctuations maintain a non-zero c-electron population. 

For G <C t the quantum fluctuations produced by the hybridization dominate, funda- 
mentally altering the behaviour of the system. Discarding the Ising interaction which is 
irrelevant for most of this regime, we arrive at the same trivial model of pseudospins track- 
ing longitudinal and transverse fields as in the G ^> t LRO phase. Here, however, the 
transverse field is generally much larger than the longitudinal field: the pseudospins will 
be mostly ordered along the x-axis with (rf) ~ |sgn(/ij). In terms of the /-occupation, 
this implies that (fjCj) ^ 0, with intermediate occupation of each /-orbital < (n?) « ~n. 
The longitudinal field produces only small deviations from this homogeneous distribution of 
electron density. This is clearly a mixed- valence state (MVS). 

Our pseudospin model allows us to connect these two extreme limits for the first time: 
for V = 0.05t we present the QFKM's phase diagram as the inset in Fig. Close to 
half-filling the Ising interaction is never dominant so the model of r-pseudospins tracking 
the longitudinal and transverse fields is valid for any G. In Fig. ([TJ we define an MVS- 
LRO boundary by equating the magnitude of the transverse and longitudinal fields. For 
small n, however, the longitudinal field is negligible; discarding this term from 7i cS we find 
a transverse field Ising model (TFIM) describes the /-orbital occupation. Although the 
TFIM has an order-disorder transition Q], the shifting of the /-level lifts the system from 
criticality and so there is a crossover regime (CR) between the SEG and MVS phases. 

The details of the CR physics are controlled by the site-dependence of the transverse 
field Eq. (jlOj) . Near the SEG phase we can use the same values for S a (xj) and C a (xj) as in 
the CP analysis; for MVS behaviour both these objects vanish. We thus assume a periodic 
transverse field within the CR. In general, the periodicity is incommensurate with the lattice: 
numerical studies 17j reveal this to be qualitatively identical to a random variation of the 

n 

transverse field. The random TFIM has been studied in depth by Fisher |Tja] using a real- 
space renormalization group treatment. To relate Fisher's results to the CR, we assume 
an effective Hamiltonian 7i CR = — J J2j r / r /+i + hj T j where the values of h°- are drawn 



randomly from the cosine distribution p{h)dh = (Cn) x y 1 — (h/C) 2 , C = 4Vy Aa/a. For 
the CR, A and a are as in the SEG phase; we therefore use the values from the classical 
(V = 0) VT problem. 
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FIG. 1: Ground-state phase diagram for the symmetric FKM below half-filling for the CP (dashed 
line) and VT problem (solid line). The numerical data points (circles for VT, plus sign for CP) are 
due to Gajek et al. P|. Inset: the phase diagram for the QFKM. Phases are as described in text. 



Lowering G from the SEG phase, the CR is reached at J c \ = max \hj\: for J < J c \ the 
Ising coupling is not greater than the magnitude of the transverse field everywhere on the 
lattice. The rare regions where > J breaks the SEG ground state up into randomly 
distributed large clusters of zero and full /-occupancy. The system nevertheless retains the 
character of the SEG phase so long as J > J& = \hj\, the average value of \W\ across the 



lattice. This is reflected by the mean / — / correlations (njnj +a .) which decay exponentially 
to (n f ) 2 > (nl) 2 with correlation length £ ~ exp{— 2[ln( l 7 c2 / l 7)] 2 }. 

Decreasing J below JT C 2, the character of the CR reverses, as here J < \hj\ over most of 
the lattice. Of the SEG state, only rare clusters remain, embedded in a valence-fluctuating 
background. This MVS-like region of the CR is characterized by a different functional 
form for the correlation length £ ~ [ln(j7" C 2/ 'J')}~ 2 ', the correlation functions decay expo- 
nentially to {nl) 2 . The CR persists until J becomes comparable to the longitudinal field 
in Eq. Q: the system is then best described by the small-G MVS. In Fig. ((TJ) the up- 
per boundary of the CR is derived from the condition J c \ = max|/tj| and is given by 



G c i = 2[{Vvpir / a5 a (a)} ^ A/aa] 1 ^ 2 ; the lower boundary is the critical line G c as in the clas- 
sical VT problem. We have extended the curve G c i(n) to G > t using the VT result for 



S 



a. 

The distinguishing feature of the CR are the randomly distributed clusters of integral 
/-valence. These clusters occur where the local c-electron density departs from the MVS 
average, determining the occupation of the underlying /-orbitals via the forward-scattering 
Coulomb interaction [i.e. Eq. (jSJ)]. This coupling of the c- and /-electron densities is recog- 
nizable as a Toyozawa electronic polaron Electronic polarons in the QFKM were studied 
by Liu, who proposed that they appear intermediate between mixed- and integral-valence 
phases |]|. Since bosonization treats the forward-scattering exactly J3], our calculation rig- 
orously confirms Liu's scenario. 

In summary, we have presented the results of a novel theoretical study of the ID FKM. 
Using a non-perturbative approach, we have uncovered the physical mechanisms responsible 
for the known /-electron ordering. We derive an effective Hamiltonian for the occupancy 
of the /-orbitals which predicts the SEG-LRO transition. We also study the QFKM for 
small hybridization V <C t. For the first time we can accurately interpolate between the 
well-known CP and MVS limits: at intermediate coupling we find a crossover regime where 
electronic polaron effects are of importance. An extensive discussion of this system is to 
follow in a later paper 
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